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Abstract. In response to the comment made by Dr. Shirakura et al ( |cond- 



mat/0011235| ) , we explain that their scaling forms of the order parameter distribution 



are inadequate. We then present an appropriate scaling form of the order parameter 
distribution, which gives a good scaling plot of the order parameter distribution itself 
and also gives consistent results with our pevious results[l]. 

In a recent paper [1], we dicussed the spin-glass phase transition in the two- 
dimensional ±J Ising model. We analysed the Binder parameter, gi, and the spin- 
glass susceptibility, xsg; by finite-size scaling including the corrections to scaling. We 
concluded that the estimation of Tc is strongly affected by the corrections to scaling, so 
that there still remains the possibility that Tc = 0. In the comment by Shirakura et al[2], 
they have performed the scaling analysis of the spin-glass order parameter distribution, 
Phio)- They insisted that two scaling forms of Pl{(i) reproduce the scaling forms of ql 
and XSG in [1]; but do not give good scaling plots of the order parameter distribution 
itself. Thus they concluded that there is no possibility that Tc = 0. 

In this reply, we first explain that the scaling forms of Pl{q) proposed in [2] are 
inadequate. Then, we present an appropriate scahng form of Pl^q), which gives a good 
scaling plot of Pl{<i) itself, and gives consistent results with our previous results[l]. 

First, the integral of rhs of equation (l)(equation (2)) in [2] gives 1/(1 + c/L^) 
((1 + d/L^yf'^). Namely, both the scaling forms of Pl{q) proposed in [2] are not 
normalized correctly. This should not be ignored, since the difference from 1 is 0(1/ L'^), 
which is the same order with the difference of the ratios of peak heights for various L in 
the scaling plots of Pl{(i)- Furthermore, when the scaling forms of Phio) (equations (1) 
and (2)) in [2] are correctly normalized, we obtain that a = and h = 2c = —d, where 
a and b are the coefficients of the correction terms in the scaling forms oi ql and xsg 
(equations (16) and (17)) in [1]. Namely, the scaling forms of Phil) proposed in [2] give 
no correction term of 0(1/ L'^) in the scaling form of ql for any values of c and d, when 
they are appropriately normalized. Thus, we conclude that the scaling forms proposed 
in [2] are inadequate. 
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We did not consider the scaling form of Pl{(i) in [1], since our method only gives 
[< q"" >l,t]p- In principle, however, the correction term of [< >l,t]p for any n comes 
from the correction term of Pl{(i)- Thus, from now on, we discuss the scaling form of 
Pl{(i) including the corrections to scahng at T = Tc. 

We think that an appropriate scaling form of Pl{(i) T = can be written as 
follows: 

P,{q) = L^I-P{qL^I-){l + (1) 

where P{x) and f{x) are some scaling functions. Pl{(i) is normalized when we take k 
as 

j P{x)f{x)dx. (2) 

Now, wc show that a simple form of f{x) gives a good scaling plot of Pl{(1) itself, and 
gives consistent results with our previous results[l]. We take /(a;) as follows: 

f 11 - 10 I a; I for I ,T |< 1.1 
[ for 1.1 <| a; I . 

First, we show the scaling plot of Pl{(i) with uo = 0.5, k = 1.252 and rj = 0.2[4] 
in figure 1, where we can see that the data with different L fit very well on one scaling 
function. (Here, we have used the data of Pl{(i) used in [2] and [3].) This P{x) gives 
g{T = 0) = 0.938 and a = —0.319, which are consistent with our previous results (figure 
7 in [1]). However, it gives Xsg{T = 0) = 1.07 and b — —0.729, which are different from 
our previous results (figure 8 in [1])[5]. Here, we show a scahng plot of xsG with u;' — 0.5, 
b = -0.729, u = 2.6, rj = 0.2 and Tc = in figure 2. We find that the data fit rather 
well on one scaling function, though the data with small linear sizes deviate from the 
scaling function. Further, Xsg(^ = 0) = 1.07 is also consistent with that in figure 2. 
Thus, we conclude that the value of b in equation (17) in [1] should be modified. It is 
noted that this modification does not change the main result in [1] that the estimation 
of Tc is strongly affected by the corrections to scaling, so that there still remains the 
possibility that Tc = 0. 

In conclusion, we have shown that the scaling forms of the order parameter 
distribution, Pl{q)-i proposed in [2] are inadequate. Further, we propose an appropriate 
scahng form of Pl{(i)- A simple example of f{x) gives a nice scaling plot of Pziq), and 
also gives consistent results with those in [1], though we have had to modify the value 
of the coefficient, b, in the scaling form of the spin-glass susceptibility, xsG- 

Acknowledgment s 

The authors would thank Dr. Shirakura for sending the detailed data in [2] and [3]. 
References 



[1] Kitatani H and Sinada A 2000 J. Phys. A 33 3545 



Comment 



3 



[2] Shirakura T Matsubara F and Shiomi M (preceding comment) cond-mat/0011235 



[3] Shiomi M Matsubara F and Shirakura T 2000 J. Phys. Soc. Jpn. 69 2798 
[4] Bhatt R N and Young A P 1988 Phys. Rev. B 37 5606 

[5] It is noted that slightly different boundary conditions have been used in [1] and [3]. Therefore, 
Pl{<i) may be different in both studies, which may lead to different values of a, b and lo. 
However, the values of g{T = 0) and xsg{T = 0) should be free from the boundary condition. 




Figure 1. A scaling plot of the order parameter distribution, Pl{<i): with uj — 0.5, 
k — 1.252 and -q = 0.2, using the scaling form of equation (1). 
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Figure 2. A scaling plot of xsg with w' = 0.5, b = -0.729, u = 2.6, t] = 0.2 and 
Tc = 0, using the scaling form of equation (17) in [1]. 



